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f*^ , We summarize the current status on constraining the density dependence of the 

symmetry energy from terrestrial laboratory measurements and astrophysical 
observations. While the value Ssym(po) and density slope L of the symmetry 
energy at saturation density po can vary largely depending on the data or 
^ , methods, all the existing constraints are essentially consistent with i?sym(po) = 

\l ■ 31 ± 2 MeV and L = 50 ± 20 MeV. The determination of the supra-saturation 

density behavior of the symmetry energy remains a big challenge. 

(N; 
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" 1. Introduction 

In nuclear physics and astrophysics, there is currently of great interest to 
_ determine the density dependence of the nuclear symmetry energy i?sym(p) 

■ that essentially characterizes the isospin dependent part of the equation of 

, state (EOS) of asymmetric nuclear matter. The exact knowledge on the 

symmetry energy is important for understanding not only many problems 
in nuclear physics, such as the structure of radioactive nuclei, the reaction 
dynamics induced by rare isotopes, the liquid-gas phase transition in asym- 
metric nuclear matter, and the isospin evolution of QCD phase diagram at 
finite baryon chemical potential, but also many critical issues such as the 
properties of neutron stars and supernova explosion mechanism in astro- 
physics.^"^ The symmetry energy may also be relevant to some interesting 
issues regarding possible new physics beyond the standard model. During 
the last decade, although significant progress has been made both exper- 
imentally and theoretically on constraining the density dependence of the 
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symmetry energy, ^'^ large uncertainties on i?sym(p) still exist, especially its 
super-normal density behavior remains elusive and largely controversial. ^^^^ 
To reduce the uncertainties of the constraints on i?sym(p) thus provides a 
strong motivation for studying isospin nuclear physics in radioactive nuclei 
at the new/planning rare isotope beam facilities around the world, such 
as CSR/Lanzhou and BRIF-II/Beijing in China, RIBF/RIKEN in Japan, 
SPIRAL2/GANIL in France, FAIR/GSI in Germany, FRIB/NSCL in USA, 
SPES/LNL in Italy, and KoRIA in Korea. 

In the present talk, we summarize the current status on constraining 
the density dependence of the symmetry energy from terrestrial laboratory 
measurements and astrophysical observations, including nuclear reactions, 
nuclear structures, and the properties of neutron stars. 

2. The symmetry energy 

The EOS of isospin asymmetric nuclear matter, given by its binding energy 
per nucleon, can be expanded to 2nd-order in isospin asymmetry S as 

E{p,S) = Eo{p) + E,y^,{p)S^ + 0{S^), (1) 

where p = pn + pp is the baryon density with pn and pp denoting the 
neutron and proton densities, respectively; 6 = (pn — Pp)/iPp + Pn) is the 
isospin asymmetry; Eo{p) = E{p,S = 0) is the binding energy per nucleon 
in symmetric nuclear matter, and the symmetry energy is expressed as 

1 d^E{p,S) ^ 

EsynAP) - ^ -Qp k=0- {^j 

Neglecting the contribution from higher-order terms in Eq. (1) leads to 
the well-known empirical parabolic law for the EOS of asymmetric nu- 
clear matter, which has been verified by all many-body theories to date, at 
least for densities up to moderate values. ^'^'^ As a good approximation, the 
density-dependent symmetry energy i?sym(p) can thus be extracted from 
the parabolic approximation as 

Esyraip) ~ E{p, 5=1)- E{p, 5 = 0). (3) 

Around the saturation density po, the nuclear symmetry energy i?sym(p) 
can be expanded, e.g., up to 2nd-ordcr in density, as 

E^y^ip) = E^yM + Lx+ + 0{x'), (4) 

where x = is a dimensionless variable characterizing the deviations of 
the density from po, and L and i^Tsym are the slope parameter and curvature 
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parameter, respectively, i.e. 
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3. The symmetry energy around the saturation density 

During the last decade, a number of experimental probes have been pro- 
posed to constrain the density dependence of the symmetry energy. Most 
of them are for the symmetry energy around the saturation density while 
a few of probes arc for the supra-saturation density behaviors. In this sec- 
tion, we summarize the present status on constraining the symmetry energy 
around the saturation density, mainly, the parameters i?syin(po) a-nd L, from 
nuclear reactions, nuclear structures, and the properties of neutron stars. 

3.1. Nuclear reactions 

Nuclear reactions, mainly including heavy ion collisions and nucleon-nucleus 
scattering, provide an important tool to explore the density dependence of 
the symmetry energy. 

3.1.1. Heavy ion collisions 

One important progress on constraining the density dependence of the sym- 
metry energy is from the isospin dependent Boltzmann-Uehling-Uhlenbeck 
(IBUU04) transport model analysis^^ on the isospin diffusion data from 
NSCL-MSU.^^ It is found that the degree of isospin diffusion in heavy- 
ion collisions is affected by both the stiffness of the nuclear symmetry 
energy and the momentum dependence of the nucleon potential. Using a 
momentum dependence derived from the Gogny effective interaction and 
the corresponding isospin dependent in-medium nucleon-nucleon scattering 
cross sections, the experimental data from NSCL-MSU on isospin diffusion 
leads to a constraint of L = 86 ± 25 MeV with Esy„,{po) =30.5 MeV,^^'^'^ 
which is shown as a solid square with error bar with a label "Iso. Diff. 
(IBUU04,2005)" in Fig. 1. It should be mentioned that the constraint 
in the original publication^''' is L = 88 ± 25 MeV and Egymipo) = 31.6 
McV, due to the application of the parabolic approximation Eq. (3) for the 
symmetry energy. This constraint is significantly softer than the prediction 
by transport model simulation with momentum-independent interaction^^ 
and in agreement with microscopic theoretical calculations. 

The isoscaling of the fragment yields in heavy ion collisions has been 
shown to be a good probe of the symmetry energy. "'^^ By analyzing 
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Fig. 1. (Color online) Constraints on -Bsym(po) and L from different experiments or 
methods. See text for details. 

the isoscaling of the fragment yields in Ar+Fe/Ca+Ni, Fe+Fe/Ni+Ni, 
Ar+Ni/Ca+Ni, and Fe+Ni/Ni+Ni reactions at Fermi energy region, within 
the antisymmetrized molecular dynamic (AMD) model, a constraint of 
£^sym(Po) = 31.6 McV and L 65 ± 25 McV has been obtained in Ref.,^* 
which is denoted by a open star with a label "Isoscaling (2007)" in Fig. 1. 

The double ratios of neutron and proton energy spectra in heavy ion col- 
lisions provide a good probe of the symmetry energy. An improved quantum 
molecular dynamics (ImQMD) transport model analysis^^ of the isospin 
diffusion data from two different observables and the ratios of neutron and 
proton spectra in collisions at E/A= 50 MeV involving -'^-'^^Sn and "'^^''^Sn 
nuclei has led to a constraint on i?sym(po) and L at 95% confidence level, 
corresponding to 2 standard deviations from the minimum x^- This con- 
straint is denoted by the region between two dotted lines with a label "Iso. 
Diff. &: double n/p (ImQMD, 2009)" in Fig. 1. A more recent ImQMD 
model analysis^^ of the isospin diffusion data from heavy ion collisions at 
lower incident energy (E/A= 35 MeV) involving "^Sn and ^^'^Sn nuclei has 
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led to a constraint of i?sym(po) = 30.1 McV and L = 51.5 MeV, which is 
shown by sohd up-triangle with a label "Iso. Diff. (2010)" in Fig. 1. 

The isospin effects of fragment transverse flows in heavy ion collisions 
are useful for extracting information on the symmetry energy. In a recent 
work,^^ the transverse flow of intermediate mass fragments (IMFs) has been 
investigated for the 35 MeV/u ™Zn +™Zn, 64zn + 64zn, and 64Ni + ^'^m 
systems. The analysis based on the AMD model with the GEMINI code 
treatment for statistically de-excitation of the hot fragments leads to a 
constraint i?sym(po) = 30.5 MeV and L = 65 MeV, which is shown by solid 
diamond with a label "Trans. Flow (2010)" in Fig. 1. 

3.1.2. Nucleon optical potential 

Experimentally, there have accumulated a lot of data for elastic scattering of 
proton (and neutron) from different targets at different beam energies and 
(p,n) charge-exchange reactions between isobaric analog states. These data 
provide the possibility to extract information on the isospin dependence of 
the nucleon optical potential, especially the energy dependence of the nu- 
clear symmetry potential. Based on the Hugenholtz-Van Hove theorem, it 
has been shown recently^^"^^ that both i?sym(/Oo) and L can be completely 
and analytically determined by the nucleon optical potentials. Averaging all 
nuclear symmetry potentials constrained by world data available in the lit- 
erature since 1969 from nucleon-nucleus scatterings, (p,n) charge-exchange 
reactions, and single-particle energy levels of bound states, the constraint 
Esymipo) = 31.3 ± 4.5 MeV and L = 52.7 ± 22.5 MeV are simultaneously 
obtained, and this constraint is indicated by the gray band with a label 
"Opt. Pot. (2010)" in Fig. 1. 

3.2. Nuclear structures 

In recent years, more and more constraints on the symmetry energy have 
been obtained from the analyses of nuclear structure properties, such as the 
nuclear mass (ground state binding energy), the neutron skin thickness, the 
nuclear isobaric analog state energies, and pygmy dipole resonances. We 
summarize these constraints in the following. 

3.2.1. Nuclear mass 

The nuclear mass data are probably the most accurate, richest, and least 
ambiguous in the nuclear data library. The Thomas-Fermi model analy- 
sis'^ of 1654 ground state mass of nuclei with N, Z > 8 has given rise to 
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£^sym(po) = 32.65 MeV and L = 49.9 MeV, which is shown by sohd star 
with a label "TF+Nucl. Mass (1996)" in Fig. 1. 

The symmetry energy coefficients asym(^) of finite nuclei with mass 
numbers A = 20 — 250 were determined from more than 2000 precisely mea- 
sured nuclear masse. With the semicmpirical connection between asym(^) 
and the symmetry energy at reference densities, i.e., Esym{pA) ~ asym(^), 
and assuming a symmetry energy with density dependence of i?sym(p) = 
Esyni{po){p/ PoV : Liu ct ul."^^ obtained a constraint at 95% confidence level 
shown as a parallelogram with short-dotted-line sides in Fig. 1, labeled 
"Nucl. Mass(2010)". 

Within the Skyrme-Hartrce-Fock (SHF) approach, it has been shown re- 
cently that a value of EgymipA) at a subsaturation reference density pA leads 
to a positive linear correlation between -Esym(po) and L.^^ Using recently 
extracted E^ymipA = 0.1 fm'^) « asy^iA = 208) = 20.22 - 24.74 MeV at 
95% confidence level from more than 2000 measured nuclear masses, Chen^^ 
obtained a constraint denoted by the region between two thick dash-dotted 
lines with a label ''E^y^{pA = 0.1 fm-3)(2011)" in Fig. 1. 

The finite-range droplet (FRDM) model has been shown to be very suc- 
cessful to describe the nuclear ground state mass. The parameters in the 
macroscopic droplet part of the FRDM model are related to the properties 
of the equation of state. Using the new, more accurate FRDM-2011a ver- 
sion, Moller et al?^ analyzed the nuclear mass of the 2003 Atomic Mass 
Evaluation (AME2003), and obtained the constraint £'sym(po) = 32.5 ±0.5 
MeV and L = 70 ± 15 MeV shown as a square box bounded by short-dash- 
dotted lines in Fig. 1, labeled "FRDM (2012)". 

In a more recent work,^^ Lattimer and Lim used the confidence ellipse 
method for nuclear mass fitting. Based on a SHF energy- density functional 
for nuclear masses, they obtained a 95% confidence ellipse for the i?syin(po)- 
L constraints shown by the thick dashed lines in Fig. 1, labeled "Nucl. 
Mass(2012)". 



3.2.2. Neutron skin thickness 

Theoretically, it has been established"'^^''^^ that the neutron skin thickness 
of heavy nuclei, given by the difference of their neutron and proton root- 
mean-squared radii, provides a good probe of i?sym(p)- The droplet model 
analyses'^"'^ on the neutron skin sizes measured in 26 antiprotonic atoms 
along the mass table leads to the constraint i5sym(po) = 28 — 35 MeV and 
L = 30 — 80 MeV shown as a square box bounded by dash-dotted lines in 
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Fig. 1, labeled "DM+N-Skin (2009)". 

A microscopic SHF analysis'^^ on the neutron skin thickness of Sn iso- 
topes has led to a set of constraints corresponding to 95% confidence levels, 
shown as a region bounded by two short-dashed curves in Fig. 1, labeled 
"SHF+N-Skin (2010)". 

3.2.3. Nuclear isobaric analog state energies 

The nuclear isobaric analog state (IAS) energies are believed to provide a 
particularly clean and useful probe of the symmetry energy since the am- 
biguities in the determination of the symmetry energy of finite nuclei from 
binding energies caused by the Coulomb term can be removed.''^ By fit- 
ting the available data on the IAS and using the droplet surface symmetry 
energy, Danielewicz and Lee'^'* obtained the constraint shown as a parallel- 
ogram bounded by dashed lines in Fig. 1, labeled "lAS+LDM (2009)". 

3.2.4. Pygmy dipole resonance 

The experimentally observed pygmy dipole (El) strength might play an 
equivalent role as the neutron rms radius in constraining the symme- 
try energy. "^^ Excess neutrons forming the skin give rise to pygmy dipole 
transitions at excitation energies below the giant dipole resonance, and 
such transitions could represent a collective vibration of excess neutrons 
against an isospin symmetric core. Comparing the measured pygmy dipole 
strength in i30,i32gj^ ^j^g^^ obtained within a relativistic mean-field ap- 
proach, Klimkiewicz et al.^^ obtained the constraint -Esym(po) = 30.2 — 33.8 
MeV and L = 28.1 — 58.1 MeV shown as a square box bounded by 
think sohd lines in Fig. 1, labeled "PDR (2007)". Another analysis^^ 
on the measured pygmy dipole strength in ^®Ni and ^'^^Sn within the rel- 
ativistic and non-relativistic mean-field approaches leads to the constraint 
Esymipo) = 31.0 - 33.6 MeV and L = 49.1 - 80.5 MeV shown as a square 
box bounded by dash-dot-dottcd hues in Fig. 1, labeled "PDR (2010)". 

3.3. The properties of neutron stars 

Astrophysical observations of neutron star masses and radii provide impor- 
tant probe for the equation of state of neutron-rich matter. In particular, 
neutron star radii are strongly correlated with neutron matter pressures 
around the saturation density.^® 

In a recent work,'^^ Steiner and Gandolfi demonstrated that currently 
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available neutron star mass and radius measurements provide a signifi- 
cant constraint on the EOS of neutron matter. Using a phenomenological 
parametrization for EOS of neutron matter near and above the saturation 
density with partial parameters determined by the quantum Monte Carlo 
calculations, they obtained a constraint of i?sym(po) = 31.2 — 34.3 MeV and 
L = 36 — 55 McV at 95% confidence level based on Bayesian analysis, -^^'^^ 
and this constraint is shown as a square box bounded by thick solid lines in 
Fig. 1, labeled "NStarl (2012)". More recently, Lattimer and Lim^^ per- 
formed a similar Bayesian analysis of the available neutron star mass and 
radius measurements, they obtained a constraint of -Bsyni(po) = 28 — 34 
MeV and L = 40 — 60 MeV shown as a square box bounded by thick 
short-dahsed lines in Fig. 1, labeled "NStar2 (2012)". 

Besides neutron star mass and radius, other properties of neutron stars 
may also put constraints on the symmetry energy. For example, the binding 
energy of neutron stars, "^^ the frequencies of torsional crustal vibrations'*^''''^ 
and the r-model instability window^^ all consistently favor L values less 
than about 70 MeV. 

3.4. Discussions 

In Fig. 1, we include totally 18 constraints on L and i?sym(po) described 
above. Obviously, it cannot be that all the constraints are equivalently re- 
liable since some constraints do not have overlap. It should be stressed 
that the symmetry energy cannot be measured directly and each constraint 
shown in Fig. 1 is based on a certain theoretical model with some approxi- 
mations or special assumptions. 

We would like to highlight two constraints in Fig. 1, i.e., "SHF+N-Skin 
(2010)" and "i?sym(pA = 0.1 fm~'^)(2011)" since both constraints are 
based on the same SHE analysis with 95% confidence. The two constraints 
are re-plotted in Fig. 2. It is very interesting to see that while the con- 
straint "i?sym(pA = 0.1 fm^'^)(2011)" indicates a linear positive correla- 
tion between L and i?sym(po), the constraint "SHF+N-Skin (2010)" dis- 
plays a negative correlation. Actually, only the constraint "SHF+N-Skin 
(2010)" among the 18 constraints displays such negative correlation. This 
interesting feature makes the constraint "SHF+N-Skin (2010)" particu- 
larly important as combing it with other constraints will significantly im- 
prove the constraint on L and i?sym(po)- It is interesting to see that the over- 
lap of "SHF+N-Skin (2010)" and ''E^y^{pA = 0.1 fm-3)(2011)" is con- 
sistent with all the other constraints shown in Fig. 1 except "lAS+LDM 
(2009)". The latter neglected the higher-order density curvature contribu- 
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Fig. 2. (Color online) Contour curves in the -Bsym(po) -L plane for -Esym{p = 0.1 
fm~^) from SHF calculations. The region between the two thick solid lines represents 
the constraint obtained with 20.22 MeV < Esym(pA = 0.10 fm-^) <24.74 MeV while 
the region between the two thick dashed lines is the constraint from the SHF analysis of 
neutron skin data of Sn isotopes within a 2a uncertainty.^^ The shaded region represents 
the overlap of the two constraints. Taken from Ref.'^^ 

tion of the symmetry energy and its inclusion may reduce the L value^^ 
(See also Ref.^^). 

The positive correlation between L and £'sym(po) from "i?sym(/OA = 0-1 
fm~'^)(2011)" has been clearly demonstrated in Ref.,^^ This feature im- 
plies that nuclear mass fitting should lead to positive correlation between 
L and £'sym(/3o) since £'sym(/3A = 0.1 fm~^) reflects the symmetry energy 
of finite nuclei, which is demonstrated by the nice agreement between the 
constraints ''E^y^{pA = 0.1 fm-3)(2011)" and "Nucl. Mass(2012)" . 

The negative correlation between L and i?sym(Po) from "SHF+N-Skin 
(2010)" can be understood from the fact that the neutron skin thickness is 
determined by the neutron and proton pressure difference at sub-saturation 
density, namely, the density slope of symmetry energy at sub-saturation 
density (rather than L), which increases with both L and i5sym(po)-^^ 

From Fig. 1 and Fig. 2, we can see that while £'sym(/3o) and L can vary 
largely depending on the data or methods, all the constraints are essentially 
consistent with -Bsym(/Oo) = 31 ± 2 MeV and L = 50 ± 20 MeV. 
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4. The symmetry energy at supra-saturation densities 

While significant progress has been made on constraining the symmetry en- 
ergy around the saturation density, the supra-saturation density behavior of 
the symmetry energy remains elusive and largely controversial. FOPI data 
on the 7r^/7r+ ratio in central heavy- ion collisions at SIS/GSI energies favor 
a quite soft symmetry energy at p > 2pQ from the isospin and momentum 
dependent IBUU04 model analysis^ while an opposite conclusion has been 
obtained from the improved isospin dependent quantum molecular dynam- 
ics (ImlQMD) model analysis. It should be mentioned that the ImlQMD 
model analysis did not consider the energy dependent symmetry potential 
and it cannot explain qualitatively the isospin fractionation phenomenon 
observed in heavy ion collisions. "'^'^^ A further careful check is definitely 
needed to understand the model dependence. 

In a more recent work, Russotto et al}^ analyzed the elliptic-flow ratio of 
neutrons with respect to protons or light complex particles from the existing 
FOPI/LAND data for ^^^Au -I- ^^'^Au collisions at 400 MeV/nucleon within 
the UrQMD model, and they obtained a moderately soft symmetry energy 
with a density dependence of the potential term proportional to {p/ Po)'^ 
with 7 0.9 ± 0.4. 

Besides using heavy ion collisions to constrain the supra-saturation den- 
sity behavior of the symmetry energy, it has been proposed recently^* that 
the three bulk characteristic parameters ii^sym(po), L and Ksym essentially 
determine the symmetry energy with the density up to about 2po- This 
opens a new window to constrain the supra-saturation density behavior of 
the symmetry energy from its density behaviors at the saturation density. 

5. Summary 

Significant progress has been made both experimentally and theoretically 
on constraining the density dependence of the symmetry energy after 
more than one decade of studies in the community Although the values 
of -Esym(po) and L can vary largely depending on the data or methods, 
all the constraints obtained so far from nuclear reactions, nuclear struc- 
tures, and the properties of neutron stars are essentially consistent with 
Esym{po) = 31 ± 2 MeV and L = 50 ± 20 MeV. More high quality data and 
more accurate theoretical methods are needed to further reduce the theoret- 
ical and experimental uncertainties of the constraints on i?sym(p) around the 
saturation density. In contrast, the determination of the supra-saturation 
density behavior of the symmetry energy is still largely controversial and 
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remains a big challenge in the community. 
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